Orbital order in the low-dimensional quantum spin system TiOCl probed by ESR 
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We present electron spin resonance data of Ti 3+ (3d 1 ) ions in single crystals of the novel layered 
quantum spin magnet TiOCl. The analysis of the g tensor yields direct evidence that the d xy orbital 
from the ti g set is predominantly occupied and owing to the occurrence of orbital order a linear spin 
chain forms along the crystallographic b axis. This result corroborates recent theoretical LDA+U 
calculations of the band structure. The temperature dependence of the parameters of the resonance 
signal suggests a strong coupling between spin and lattice degrees of freedom and gives evidence for 
a transition to a nonmagnetic ground state at 67 K. 

PACS numbers: 71.27.+a, 76.30.Fc, 75.10.Jm 



Transition metal (TM) oxides with low-dimensional 
structural elements provide a fascinating 'playground' to 
study novel phenomena such as high-temperature super- 
conductivity, spin-charge separation, spin-gap states and 
quantum disorderiS*^. Until recently, the emphasis has 
been put on Cu-based oxides, where a Cu 2+ (3d 9 ) ion has 
a single hole in the e g orbitals with spin 5 = 1/2, and its 
orbital momentum is almost completely quenched by the 
crystal field. The ions at the beginning of the TM ele- 
ments row, like Ti 3+ and V 4+ , have, in contrast, a single 
d-electron which occupies one of the tig orbitals. Because 
these orbitals are much less Jahn- Teller active, their near 
degeneracy may yield more complicated physics, involv- 
ing not only the spin and charge, but also the orbital 
sector—. As an example, the three-dimensional cubic per- 
ovskite LaTiC>3 has been proposed to realize a quan- 
tum orbital liquid 5 ' 6 . However recent x-ray and neu- 
tron structural data suggest the ordering of the orbitals 7 . 
The structural dimensionality is reduced in TiOCl, where 
[TiO^Cy octahedra are arranged in bilayers separated 
from each other along the c axis (Fig. [IJi). In fact, for 
quite a while this compound has been considered as a 
2D antiferromagnet, an electron analog to the high-T c 
cuprates 9 , owing to an almost T-independent magnetic 
susceptibility reported in Ref. H. However, very recently 
TiOCl has emerged in an entirely new light as a ID 
antiferromagnet 19 and is proposed as the second example 
of an inorganic spin-Peierls compound after CuGeO^ 11 . 
LDA+U band-structure calculations^ , suggest ordering 
of the t2 g orbitals in TiOCl which produces quasi-lD an- 
tiferromagnetic (AF) S = 1/2 chains. This calculation 
favors the occupancy of the d xy orbitals (Fig. |TJ)) which 
form a uniform chain along the b axis. A transition to 
a non-magnetic state at T c = 67 K has been observed in 
the static magnetization 10 . Remarkably, NMR data re- 
veal the pre-existing pseudo spin-gap already above T c 
which is ascribed to strong orbital fluctuations^. 

In this paper we present electron spin resonance (ESR) 
data of Ti 3+ (3d 1 ) in single crystals of TiOCl. By mea- 
suring the anisotropy of the g factor and comparing it 
with our theoretical estimates in the framework of the 
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FIG. 1: (a) Crystal structure of TiOCl; (b) sketch of the t 2g 
orbitals (a, f3 and 7) in the [TiCUCy octahedron. The local 
coordinate frame {xyz} is chosen so that z|| a axis, and the x 
and y axes are rotated by 45° with respect to the c and b axes, 
respectively. Note that in the LDA+U calculation of Seidel et 
al*— one- dimensional bands are formed from overlapping xy, 
yz, and xz orbitals, respectively, whereas on a single cluster 
level the a = xy, (3=(xz+yz) /v2 and 7= (xz— yz) / \/2 states 
are the eigenstate orbitals, respectively. 



angular overlap model we conclude that the single d elec- 
tron occupies the d xy orbital which lies in the be plane. 
This result suggests a formation of a spin-1/2 chain along 
the b direction, owing to the overlap of the orbital states, 
and supports recent LDA+U calculations 10 . The ESR 
signal vanishes at T c = 67 K signalling the transition to 
a non-magnetic ground state. A pronounced dependence 
of the linewidth and the g factor on temperature suggests 
a strong coupling of spins to the lattice which may play 
an important role for the opening of the spin gap. 

Single crystals of TiOCl have been prepared from 
TiCbj (Aldrich) and Ti02 (Kronos Titan) according to 
experimental details given in Ref. Il-Sl The purity of the 
product was checked by x-ray powder diffraction at 293 K 
and 10 K. The latter measurement has been carried out 
in order to check for a possible structural phase transi- 
tion, which we do not fincU^. Both diffractograms could 
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FIG. 2: (a) Static susceptibility x(T) of TiOCl together with 
the ESR spin susceptibility and corresponding fits (see text); 
(b) and (c) ESR linewidth AH and g factor, respectively. 



be indexed in the orthorhombic space group Pmmn with 
a = 3.789(1) A, 6 = 3.365(1) A, c = 8.060(3)A at 293 K 
and a = 3.7946(3) A, 6 = 3.3584(2) A, c = 8.057(1) A at 
10 K. It is interesting to note that the a axis is longer 
and the 6 and c axes are shorter at lower temperatures, 
which means that the already distorted [TIO4CI2] octa- 
hedra are slightly less compressed at 10 K. The static 
magnetic susceptibility x(T) has been measured in a 
Faraday balance in a field of 1 Tesla in the tempera- 
ture range 10—550 K (Fig. Ek) ■ Within the experimental 
accuracy x(T) is isotropic and shows a behavior similar 
to that reported in Ref. 0. The fit of the data with 
the formula x(T) = xCHidaf + xC^Ocurie + Xo, where 
x(^)idaf is the susceptibility of a uniform ID Heisen- 
berg antiferromagnefei^, x(^) Curie is the Curie term, and 
Xo is the sum of the diamagnetic susceptibility of Ti- 
OCl and the Van-Vleck susceptibility of Ti 3+ , yields the 
number of free 5=1/2 spins (paramagnetic defects), re- 
sponsible for the low-T Curie upturn of x(^) °f order 
0.6 %, and xo = 2.3 x 10~ 4 emu/molei&. After the sub- 
traction of x(^) Curie and xo from the raw data one sees 
that the model of a ID antiferromagnet with the nearest 
neighbor exchange J = 676 K and the g factor of 1.91 
(see below) describes the data in the high-T regime quite 
well (Fig. [5Jl). Below 130 K the susceptibility rapidly 
decreases and approaches zero at ~ 60 ]^Mtl^ii^_ 

ESR has been measured using a Bruker spectrometer 
at X-band frequency 9.48 GHz and at temperatures be- 
tween 2 and 300 K. A single resonance line of a Lorentzian 
shape has been observed. The intensity of the ESR sig- 
nal I is proportional to the susceptibility of the resonat- 
ing spins 20 . A comparison of I of TiOCl with that of 
a reference sample^ gives evidence that within the ex- 
perimental uncertainty practically all spins contributing 
to the static susceptibility participate in ESR. The In- 
dependence of the spin susceptibility ^ ESR obtained from 
the integrated intensity I(T) is similar to that of the 
static x (Fig- I2t) 21 - I n particular, x R drops sharply 
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FIG. 3: Angular dependence of the g factor at T=100 K: (a) 
(1-10) plane; (b) (001) plane. 



almost to zero at T c = 67 K signalling the transition to a 
non-magnetic state. In addition a small kink can be seen 
at 90 — 95 K, which is ascribed in Ref. to the onset of 
the structural instability. A small signal arising at T < T c 
is obviously due to a small amount of paramagnetic im- 
purities in the samples (e.g. Ti 3+ in structural defects). 
The bulk signal observed above T c is anisotropic and ex- 
hibits a strong temperature dependence of the g factor 
and the linewidth Ail (Fig. 03 and c). 

A representative example of the dependence of g on 
the direction of the external magnetic field H is shown 
in Fig. [3] One notices that the c axis is a symmetry 
axis of the g tensor, where g reaches its minimum value 
(Fig. |3^). The g factor in the ab plane is larger and 
also anisotropic, as expected in the orthorhombic sym- 
metry (Fig. |2t>) : The situation is typical for a Ti octahe- 
dron strongly compressed along the z axis with smaller 
distortions in the equatorial plane, where the xy orbital 
ground state is realized^ 3 .. On the quantitative level we 
compare the experimental data with the results of the 
angular overlap model (AOM). The AOM has been a 
tool to rationalize spectroscopic and magnetic properties 
of TM and rare-earth metal ions in various ligand fields 
and provides a mathematical approach to calculating rel- 
ative energies of molecular orbitals of a TM complex from 
the overlap of the central-atom orbitals with the ligand 
orbitals^i. We calculated the g values for TiOCl using the 
program package CAMMACSSi. The best fit of the g 
tensor with experimental data (Table clearly identifies 
the d xy orbital as the ground state orbital. The relative 
energies of the orbital states and the bonding parame- 
ters used in the calculation are listed in Table [H] Even 
if one neglects a plausible anisotropy of the it interac- 
tions (Model A) the AOM predicts correctly a minimum 
value of the g factor along the c axis and larger values 
in the ab plane (g c < gb < g a )- The agreement becomes 
even quantitative if one accounts for the anisotropy of the 
7r bonding, which is quite reasonable owing to the fact 
the ligands are all bridging, e.g. [Ti0 4 /4Cl 2 /2] (Model 
B)Sa. In the LDA+U calculation of Seidel et alJii. the 
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TABLE I: Observed and calculated g tensor at T 


= 293 K. (A 


- isotropic, B - anisotropic n interactions). 




experiment Model A 


Model B 


g a 2.01 


1.946 


1.976 


g b 1.96 


1.935 


1.959 


g c 1.91 


1.926 


1.911 


TABLE II: AOM parameters for [Ti0 4 Cl 2 ] complex and rel- 


ative energies of the d orbital states. 
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C 2800 cm -1 






e CT 9500 cm -1 , d(Ti-O) = 1.95 A 






e w 2400 cm -1 


z 2 


20769 cm -1 


e CT 4900cm -1 , d(Ti-O) = 2.25 A 


2 2 

x -y 


14891 cm -1 


1200 cm -1 


(xz+yz)/V2 


5900 cm -1 


e CT 5700cm -1 , d(Ti-Cl) = 2.37 A 


(xz-yz)/V2 


2456 cm -1 


e n 900 cm -1 
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k x 0.6 k y , k z 0.9 







two higher-energy ID bands are derived from unmixed 
d yz and d xz orbitals, which both make an angle of ~ 45° 
with the c axis and the ab plane. Their occupation by 
the d electron can be excluded, as the symmetry of these 
states is different from the symmetry of the g tensor. 

The above discussion strongly supports the scenario 
proposed in Ref. on the basis of the band structure 
results. The occupation of the d xy states favors the oc- 
currence of AF spin chains along the b axis owing e.g. 
to the direct overlap of the d xy orbitals (Fig. The 
neighboring chains in the [TiC^Cla]^, bilayer are shifted 
in a staggered fashion by a/2 and 6/2 along the a and 
b axes, respectively (Fig. QJi), so that the inter-chain ex- 
change interaction is expected to be weak and frustrated. 
Thus almost perfectly isolated spin- 1/2 chains with a rel- 
atively strong AF exchange J ^700 K, as deduced from 
the x data, can be realized. The spin chains undergo 
the transition to a spin-gap state at T c = 67 K. How- 
ever, NMR datai£ indicate that this is not a conventional 
spin-Peierls transition as observed e.g. in CuGeO^ii be- 
cause of a much larger spin gap A gap as compared with 
T c (2A gap /fc s T c ~13) 14 . Probably additional orbital de- 
grees of freedom play an important role in TiOCl. In- 
deed, the components of the g tensor exhibit a strong 
T dependence (Fig. [2t), suggesting that the energy of 
the orbital states alters appreciably with T. This is a 
rather unusual feature; for example, in CuGeC>3 as well 
as in NaV205 which is another well-known low-D spin- 
gap TM oxide, the g values do not change with 722*21. 



A considerable decrease of the g anisotropy at low T 
signals a reduction of the distortion of the Ti complex, 
in agreement with our x-ray data. Such an apprecia- 
ble coupling of the spin to the lattice should affect the 
T dependence of the ESR linewidth as well. In a con- 
centrated paramagnet spin-spin interactions produce a 
finite ESR linewidth owing to the anisotropic part of the 
superexchange^S TL 1 = ^2 SiASj. This interaction yields 
a second moment of the line Mi (Ag/g) 2 J Here 

Ag is the deviation of g from the spin-only value of 2, and 
J is the strength of the isotropic Heisenberg exchange 
TL = J ^2 SiSj which in the 3D case narrows the signal so 
that its width reduces to AH ~M%/J Mdl If the suscep- 
tibility x(T) deviates from the Curie law, M 2 may acquire 
a T dependence proportional to XCurie/xC^ 1 ) ^ Taking 
Ag/<7~0.05 and J~680 K we obtain AH of the order of 
30 Oe. The much smaller value of this rough estimate as 
compared with the experimental data for T>T C (Fig. 03) 
may be in part due to the fact that it neglects the pecu- 
liarities of the bonding geometry in the ID spin chain 
which in certain cases may considerably boost AH? 7 
More serious is the discrepancy of the T dependence of 
AH. Because XCurie/x(^) an d the average value of Ag 
both increase with lowering T, the width AH is expected 
to increase, too. However, experimentally one finds the 
opposite behavior (Fig. 12b), suggesting that other mech- 
anisms of spin relaxation e.g. via orbital and lattice de- 
grees of freedom have to be considered^. The interplay 
between orbital and spin fluctuations in TiOCl has been 
proposed in the discussion of the NMR results^. It is ar- 
gued that the opening of the pseudo-spin gap at T* ~ 135 
K is related to the suppression of the spin-Peierls transi- 
tion caused by fluctuations of the orbital states. Remark- 
ably, in the temperature interval T c < T < T* the ESR 
linewidth levels off at a minimum value of ~ 150 — 170 
Oe before dropping down at T c . This implies that the 
spin fluctuations in this temperature regime are strongly 
suppressed as expected in the pseudo-gap regime. The 
spin dynamics recovers above T* resulting in the increase 
of AH. An additional T-dependent contribution to AH 
could arise due to the spin-phonon coupling, which may 
be significant as suggested by the strong sensitivity of the 
g factor to the change of the lattice parameters^. 

In summary, we have studied electron spin resonance of 

Xi 3+ 

ions in single crystals of the novel low-dimensional 
spin magnet TiOCl. The analysis of the g tensor justifies 
the scenario that in an apparently 2D structure uniform 
spin 5 = 1/2 chains are formed along the b direction. 
The bulk ESR signal vanishes at T c = 67 K evidencing the 
transition to a nonmagnetic, possibly a spin-Peierls state. 
The T dependence of the g values and the linewidth sug- 
gests that orbital and lattice degrees of freedom may play 
a key role in the magnetic properties of TiOCl. In par- 
ticular, strong spin and probably also orbital fluctuation 
effects above T c may be responsible for a peculiar tem- 
perature dependence of the ESR parameters. 
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